Excess current in point contacts on two-band superconductor MgB2 in magnetic field 
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A series of I{V) characteristics and bias-dependent differential resistances dV/dI{V) of the point 
contacts made of the single crystal of two-band superconductor MgB2 were measured in magnetic 
fields up to 9T. The magnetic field dependences of the excess current in the /(V) curves were ob- 
tained and analyzed using Koshelev and Golubov's [Phys. Rev. Lett. 90, 177002 (2003)] theoretical 
results for the mixed state of a dirty two-band superconductor. Introducing a simple model for the 
excess current in the point contact in the mixed state, our data can be qualitatively described us- 
ing the theoretical magnetic field dependence of the superconducting order parameter of the a and 
TT-bands and the averaged electronic density of states in MgB2. 

PACS numbers: 74.45.-(-c, 74.70Ad, 73.40.Jn 



I. INTRODUCTION 



II. SAMPLES AND MEASUREMENTS 



X 
S3 



MgB2 has attracted considerable attention due to its 
high Tc (w40K) which is due to specific multi-band elec- 
tronic structure (see Ref.i and Refs. therein). Here 
superconductivity develops in two-dimensional cylindri- 
cal sheets of the Fermi surface constituted from the a- 
band and three-dimensional tubular Fermi surface from 
the TT-band. The commonly accepted explanation for 
high-Tc superconductivity in MgB2 is connected with the 
strong interaction between the charge carriers and the 
i?2g phonon mode,^ caused by the antiparallel vibrations 
of the atoms in the honeycomb-like boron planes. The 
superconducting order parameter is distributed over the 
Fermi surface of MgB2, being Ao- « 7 meV for the cr-band 
and Att ~ 2meV for the 7r-band.^ The different prop- 
erties of the a and tt bands with different anisotropics 
result in very peculiar and interesting physical charac- 
teristics of MgB2 both in the normal and especially in 
the superconducting stateii*^ 

Point-contact (PC) spectroscopy is one of the straight- 
forward methods which has proved the existence 
the double-gap superconductivity in MgB2 (see, e.g., 
Refs.^iSi^). Although much was done to understand the 
remarkable superconducting state in this diboride, some 
issues should be clarified, such as the behavior of the su- 
perconducting gaps in the magnetic field and the rapid 
depression of the superconducting features of the 7r-band 
in the PC conductivity by a moderate field of about 1 T. 
In this study we carried out PC measurements both of the 
order parameter and the excess current versus the mag- 
netic field in MgB2 and analyzed them using the theoryS, 
for the mixed state of a dirty two-band superconductor. 
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Our experiments were carried out on the MgB2 single- 
crystal samples investigated earlier in Ref . . Their 
properties are described elsewhere. The crystals were 
small thin plates (flakes) of a sub-millimeter size. They 
were glued with silver epoxy to the sample holder at one 
of their side faces. The noble metal counter electrode was 
brought into a gentle touch with the opposite side face 
of the crystal in liquid helium. Thus we tried to make a 
contact mainly along the ah plane, for which both the a 
and TT-band superconducting features should be present 
in the PC characteristics. The magnetic field was applied 
along the c-axis. 

The I{V) characteristic and the differential resis- 
tance dV/dI{y) were recorded through sweeping DC 
current / on which small AC current i was superim- 
posed and through measuring the alternating voltage 
Vi oc dV/dI{V) using the standard lock- in technique. 



III. EXPERIMENTAL RESULTS 

First we evaluated the characteristic length scales im- 
portant for this study. The contact size for the typi- 
cal resistance of 1 ^7 is estimated by the Sharvin for- 
mula d ~ pl/R {d ~ 7nm in Ref.—) irrespectively 
of the crystallographic directions. For a dirty constric- 
tion (with a very short mean free path) the use of the 
Maxwell formula d ^ p/R results in d about 0.7 nm 
and 2.6 nm (see Ref.—) for the ab plane and c direc- 
tion, respectively. Thus, the contact size can be of the 
order of magnitude or smaller than the elastic electron 
mean free path (lab — 24 nm, Ic = 6nm, according to 
Ref.—) , which means that theoretically we can approach 
the ballistic (spectroscopic) regime. In practice the spec- 
troscopic regime is confirmed by observation of the char- 
acteristic Andreev- reflection features in dV/ dl in the su- 
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FIG. 1: Magnetic field dependence of large and small super- 
conducting energy gaps at 4.2 K (solid squares and triangles, 
left y-axis) obtained by a BTK fitting of the curves shown in 
the inset. Solid curves show magnetic field dependence of the 
maximal pair potential (right y-axis) for MgB2 for Di = D2 
(see Fig.|5J. Inset: symmetrized experimental dV/dl curves 
in magnetic field (solid lines) for the single crystal MgB2-Cu 
2.2 Q, junction along the ah plane with BTK fit (thin lines) 
from Ref.J^. Two separate sets of the gap minima around 
±2.5 and ±7 mV are clearly seen in the low fields. 



perconducting state. 

The measurements of the superconducting energy gap 
in MgB2 by means of the Andreev reflection demonstrate 
two sets of the energy-gap minima in the differential re- 
sistance dV/dl at ±2.4± 0.1 and ±7.1±0.4meV, respec- 
tively; their distribution^ corresponds to the theoretical 
predictioni^ The dV/AI curves can be well fitted by the 
known BTK^'* equations (with the suitable F parameter) 
for two conducting channels with the appropriate weight 
factors y;,SiSiL& The contribution of the large gap to the 
double-gap spectra is the largest along the ah plane 
and amounts theoretically to 0.34, while it is negligible 
in the perpendicular direction (along the c axis)pi^ 

A series of dV/dl measured in magnetic field is shown 
in Fig.n(inset) with clearly seen two sets of the gap 
minima (double-gap features). The main panel in Fig.^ 
displays the magnetic field dependence of the large and 
small gaps obtained by the routine BTK fitting of the 
curves from the inset. In spite of the rapid decrease in 
the intensity of the small gap minima with the field (a 
factor of two at 1 T, see the inset in Fig.^, the small gap 
is not suppressed by low fields (< IT), and the estimated 
critical field 6-7 T (see Fig^ is much higher than that 
statedSS in Refs.^*i&. Moreover this field corresponds to 
Bc2 determined from the onset temperature in the broad- 
ened resistive transition in single crystals^l^ Similar be- 
havior of the small gap persisting up to 5 T was found in 
Ref.— , where visual disappearance of the 7r-gap features 
in the PC characteristics was attributed to broadening of 
the spectra in the magnetic field. 

Experimental data for A in Fig.^ are in good agree- 



FIG. 2: Magnetic field dependences of the maximal pair po- 
tential A and averaged electron DOS A''(0, H) at zero energy 
for two ratios of the diffusion constants Di = for the a 

(1) and TT (2) bands. By courtesy of Koshelev and Golubovi- 

ment with the magnetic field dependence of the maxi- 
mal pair potential A for MgB2 according to Ref. ^ with 
Di ~ D2 relation. Under this condition, the smaller A2 
(see Fig. 121 dashed curve) has a negative curvature in low 
and near to the critical fields, which is similar to the 
experimental A5 (see Fig.QJ. 

The Andreev refiection of electrons from the N — S 
boundary in a — c — 5 point contact (A^ stands for the 
normal metal, c is the constriction, and S is the super- 
conductor) produces an excess current lexc in the super- 
conducting state. As a result, I{V) can be writtenM at 
ey > A as I{V) = V/Rn + /cxc, where 

/e.c « A/ei?^v (1) 

with A being the superconducting energy gap and Rn 
being the PC resistance in the normal state. 

We determined the reduced Iexc{B) / Iexc{B = 0) either 
from the raw I{V) or mostly by integrating the reduced 
R-^dV/dI{V) or RdI/dV{V) after subtracting the nor- 
mal state background. The use of the I{V) curves pro- 
vides less accurate results because of the resistance in- 
stability for some contacts during the field sweep. 

For all investigated contacts the excess current depends 
on the magnetic field with the overall positive curvature 
(Fig. 01 . At first, Iexc{B) decreases abruptly and then 
more slowly above 1 T. This corresponds to the afore- 
mentioned drastic suppression of the dV/dI{V) small-gap 
minima intensity by a magnetic field about 1 T and to the 
persistence of the residual superconducting minima on a 
further increase in the magnetic field. This is quite dif- 
ferent from what is expected for lexc according to 
taking into account that A(i3) has, as a rule, a negative 
curvature (see FigEJ. 

Generally, the excess current, just like the net current 
through the contact, is proportional to the contact area 
A (see, e. g., Eqs. (18), (19) in Ref. "), i. e. lexc ^ Ax A. 
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Assuming that lexc is caused by the superconducting 
component of the current, which in the mixed state is 
determined by the superconducting volume (part) in the 
PC area, the coefficient A in the above-mentioned for- 
mula should be effectively reduced by the contribution of 
the nonsupcrconducting volume. Koshelev and Golubo\«^ 
calculated the magnetic field dependence of the electron 
density of states (DOS) N{E) at E=0 averaged over the 
vortex unit cell. We assume that in the first approxima- 
tion the averaged iV(0, B) can be regarded as a measure of 
the nonsupcrconducting region in the mixed state, and, 
vice versa, (1 — A^(0, B)) characterizes the superconduct- 
ing phase portion. Thus, according to above reasoning, 
we suppose that 

lexc oc WLil ~ Nl{0, B))Al + ws{l - Ns{0, B))As (2) 

with the weight factors wl and ws for large and small 
gaps, respectively {wl + ws = 1), Nl,s{0,B) and 
^{B)l,s are taken from Fig.|21 where L = 1, S = 2. 
Here, I^xc is presented^^ by iPJ as a sum of two contribu- 
tion from the a and tt bands weighted by wl,s factors, 
which along with A were derived from the fitting of the 
Andreev-reflection features in dV/dl at zero field. 

Figure|21 shows the magnetic field dependence of lexc 
for three different types of dV/dl with visually distinct 
small, large and double gap structures. Here, we also 
plotted lexc calculated by lj2Jl. The most unexpected 
result is the observation of dV/dl (see Fig.l^Jb)) with 
large-gap minima only, which is discussed in Ref. . The 
large-gap minima are only slightly affected by magnetic 
fields of a few tesla. This is in contrast to the small-gap 
minima, which according to Refs.^i^ should vanish above 
1 T. The curve in the bottom panel was calculated tak- 
ing Wl from the fitting of the double-gap structure (see 
Fig.0, while in panel (a) and (b), where only the small- 
or large-gap structures are seen, wl can vary from zero to 
its highest value 0.34. — The latter wl value was used for 
the spectra in the middle panel in which only large-gap 
minima are seen. Evaluation of wl for the spectrum in 
panel (a) inset is not straightforward because the spikes 
present at higher energies mask the large gap features. 
However, selected wl =0.15 and wl = 0.34 fit reason- 
ably both zero-field dV/dl and Iexc{B) curves in panels 
(a) and (b), respectively. 

In general, the experimental and calculated curves 
agree well in spite of the simple empiric model adopted 
for lexc- The calculated curves describe qualitatively 
the experimental data, especially the abrupt initial de- 
crease in lexc, pronounced in the curves with the domi- 
nant small-gap features (Fig.|2|^a)). In Fig.O^c), showing 
two curves with different D1/D2 ratios, it is seen that 
better correspondence can be reached with an interme- 
diate ratio between Di values. The same occurs for the 
data in Fig.E^b) for the large-gap spectrum, though the 
experimental data for large-gap and small-gap spectra 
are closer to the cases with Di = D2- This fact testifies 
that in this case the a band is closer to the clean limit, 
also mentioned in Ref.—. 
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FIG. 3: lexc vs magnetic field behavior for 3 PCs between the 
MgB2 single crystal and Cu at 4.2 K: (a) junction with small- 
gap minima and normal state resistance i?=1.5 Q., (b) junction 
with large-gap minima from Ref. R—&.7Q,, (c) junction 
with double-gap structure from Fig.Q R—2.2 S7. The lines 
present calculated Iexc{B) according to l|5Jl with wl, Di and 
D2 shown in each panel. The error bars show difi'erence be- 
tween Iexc{B) calculated from the experimental (insets, solid 
curves) and BTK (insets, dashed curves) dV/dl. Insets: zero- 
field experimental (solid curves) and BTK (dashed) dV/dl 
curves for each contact with As=2.55meV and AL=7.8meV 
(a), As=2meV and AL=7.5meV (b), and As=2.55meV and 
AL=7.1meV (c). Vertical bars are drawn to show relative 
dV/dl value with respect to the contact resistance R. 



IV. CONCLUSIONS 

We have investigated the excess current in PCs based 
on single crystals of MgB2 in the magnetic field and pro- 
posed an empiric model for the excess current behavior 
in the mixed state. Using this model, a magnetic field 
dependence of the excess current was obtained, which is 
in qualitative agreement with experimental data. The 
anomalous magnetic field dependence of the excess cur- 
rent in PCs reflects the specific two-band structure of the 
superconducting order parameter and the spatial DOS 
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behavior in MgB2. Thus, Iexc{B) can be used for eval- 
uation of the diffusivity in both bands on the base of 
Koshelev and Golubov's theoryi^ 
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